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We study a number of different ingredients, related to long range order observed in lattice QCD 
simulations, using a simple "deformed QCD" model. This model is a weakly coupled gauge theory, 
which however has all the relevant crucial elements allowing us to study difficult and nontrivial 
questions which are known to be present in real strongly coupled QCD. Essentially, we want to un- 
derstand the physics of long range order in form of coherent low dimensional vacuum configurations 
observed in Monte Carlo lattice simulations. 



I. INTRODUCTION AND MOTIVATION 

The main motivation for this work is the recent Monte 
Carlo studies of pure glue gauge theory which have re- 
vealed some very unusual features. To be more specific, 
the gauge configurations display a laminar structure in 
the vacuum consisting of extended, thin, coherent, locally 
low-dimensional sheets of topological charge embedded in 
4d space, with opposite sign sheets interleaved, see orig- 
inal QCD lattice results [1-4]. A similar structure has 
been also observed in QCD by different groups [5-9] and 
also in two dimensional CP N ~ X model [10]. Furthermore, 
the studies of localization properties of Dirac eigenmodes 
have also shown evidence for the derealization of low- 
lying modes on effectively low-dimensional surfaces. The 
following is a list of the key properties of these gauge 
configurations which we wish to study: 

1) The tension of the "low dimensional objects" van- 
ishes below the critical temperature and these objects 
percolate through the vacuum, forming a kind of a vac- 
uum condensate; 

2) These "objects" do not percolate through the whole 
4d volume, but rather, lie on low dimensional surfaces 
1 < d < 4 which organize a coherent double layer struc- 
ture; 

3) The total area of the surfaces is dominated by a 
single percolating cluster of "low dimensional object"; 

4) The contribution of the percolating objects to the 
topological susceptibility has the same sign compared to 
its total value; 

5) The width of the percolating objects apparently van- 
ishes in the continuum limit; 

6) The density of well localized 4d objects (such as 
small size instantons) apparently vanishes in the contin- 
uum limit. 

It is very difficult to understand the above proper- 
ties using conventional quantum field theory analysis. 
Indeed, the QCD lattice results [1-8] imply that the 
topological density distribution is not localized in any 
finite size configurations such as instantons; rather the 
topological density is spread out on the surface of low- 
dimensional sheets. Such a structure can not be immedi- 
ately seen in gluodynamics, at least not at the semiclas- 
sical level. At the same time, these Monte Carlo results 
could be interpreted very nicely with a conjecture that 
the observed structure is identified with the extended D2 



branes in holographic description[ll-14]. 

One of the key elements of this conjecture is assump- 
tion that the tension of the D2 branes vanishes below 
the QCD phase transition T <T C such that an arbitrary 
large number of these objects can be formed. Vanishing 
tension in the dual description in the confined phase is 
a result of the Hawking-Page phase transition [15] when 
the D2 brane shrinks to the tip of a cigar type geometry. 

The second key element in identification of the struc- 
ture observed on the lattice [1-8] with the holographic 
description in terms of the D branes is the assumption 
that the topological density distribution which is orig- 
inally localized in well defined DO branes (instantons), 
somehow spreads out along extended D2 branes as a re- 
sult of strong interaction between D0-D2 branes, lead- 
ing to their binding. Such a picture was basically mo- 
tivated, as mentioned in [12-14], by the structure which 
emerges in supersymmetric field theories [16] where the 
relevant dynamics can be indeed formulated in terms of 
the strongly bound D0-D2 configurations. 

In this paper, we investigate precisely the second idea 
above in the framework of a "deformed QCD" developed 
in [17]. The deformation allows us to bring the gauge 
theory into a weakly coupled regime wherein calculations 
can be performed in theoretically controllable manner. In 
spite of the great deal of analytic control provided, the de- 
formed theory preserves many of the relevant structures 
present in strongly coupled QCD including confinement, 
degeneracy of topological sectors, and the correct non- 
trivial 6 dependence. Furthermore, it seems, there is no 
order parameter differentiating the weakly coupled de- 
formed regime from the strongly coupled regime, which 
reproduces undeformed QCD [17], so that the (gross) be- 
haviour of the two theories may be quite similar. 

In particular, the deformed theory exhibits two im- 
portant structures of note: first, the topological charge 
in this model is carried by the fractionally charged 
monopoles with topological charges Q — ±l/iV; and sec- 
ond, there are domain walls present in the system as a 
result of a generic 2n periodicity of the effective low en- 
ergy Lagrangian governing the dynamics. Given these 
ingredients, we would like to test the following two ideas 
which are apparently related to the configurations ob- 
served in the lattice simulations [1-8]: 

1) the domain walls form precisely a double layer struc- 
ture with opposite sign sheets of the topological charge 
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density interleaved; 

2) the monopoles and domain walls attract each 
other and the topological charge originally localized on 
monopoles spreads out along the domain walls. 

If the second occurs, there will be few well-localized 
finite sized sources carrying the topological charge. In- 
stead, the topological charge density will be spread over 
extended domain walls, which is precisely what has been 
observed in simulations [1-8]. 

We note that a similar picture of attraction between 
monopoles and domain walls was originally discussed in 
a cosmological context [18], see also the related papers 
[19-21] and references therein. The basic idea in [18] is 
that if physical monopoles and domain walls are present 
in the system, there will be an attractive force between 
them. Then, if these objects collide, the monopole's 
winding number (monopole charge) spreads out on the 
surface of the domain wall, and will be eventually pushed 
to the boundaries at infinity. This effect was suggested 
as a solution of the so-called "cosmological monopole's 
problem". In our context we do not have real physical 
monopoles and real physical domain walls in Minkowski 
space, but rather Euclidean monopoles and domain walls 
which must be interpreted as configurations describing 
the tunnelling processes in physical Minkowski space, see 
detail discussions of this point in [13]. Nevertheless, the 
formal structure of the problem and relevant features 
(such as attraction between the objects and spreading 
the magnetic charge over the surface) are very much the 
same. 

We also note that "deformed QCD" model has been 
successfully used to test some other nontrivial features 
of strongly coupled QCD such as emergence of non- 
dispersive contact term in the topological susceptibility 
[22] and the emergence of a topological Casimir behaviour 
in gauge theory with a gap [23]. In both cases the effects 
of interest are a result of the nontrivial topological vac- 
uum structure of this model. 

The structure of our presentation is as follows. In next 
section (II), we review the relevant parts of the model 
[17, 22] including the low energy description of the theory 
in terms of the sine-Gordon Lagrangian. In the following 
section (III), we construct the domain walls and explic- 
itly demonstrate the double layer structure apparently 
observed on the lattices. Finally, in Section IV we study 
the interaction of the domain walls and monopoles. 



II. DEFORMED QCD 

In the "deformed" Yang-Mills, developed in [17], an 
extra "center-stabilization" term is put into the the La- 
grangian in order to prevent the center symmetry break- 
ing that characterizes the QCD phase transition between 
"confined" hadronic matter and "deconfined" quark- 
gluon plasma. Thus we have a theory which remains 
confined at high temperature in a weak coupling regime, 
and for which it is claimed [17] that there does not exist 



an order parameter to differentiate the low temperature 
(non-abelian) confined regime from the high temperature 
(abelian) confined regime. We now proceed, in section 
II A, to review the relevant aspects of the theory. We 
then discuss, in section II B, the low-energy effective La- 
grangian which gives rise to the domain wall solutions 
mentioned earlier. 



A. Formulation of the theory 

We start with pure SU (N) Yang-Mills (gluodynamics) 
Wick rotated with Euclidean time compactified on the 
manifold R 3 x S 1 defined by the standard action 



We then add to it a deformation action, 



AS= [ d 3 x^PP(x)], 
Jr 3 L 



(i) 



(2) 



built out of the Wilson loop (Polyakov loop) wrapping 
the compact dimension, 



fl(x) 



i § dX4 A4(x,X4) 



(3) 



The "double-trace" deformation potential P [ft] respects 
the symmetries of the original theory and is built to sta- 
bilize the phase with unbroken center symmetry. It is 
defined by 



LAT/2J 

J2 an 



|tr [n n 



(4) 



Here [N/2\ denotes the integer part of N/2 and {a n } is 
a set of suitably large positive coefficients. 

In the undeformed theory the effective potential for 
the Wilson loop is minimized for Q an element of Zjv- 
The deformation potential (4) with sufficiently large {a n } 
however changes the effective potential for the Wilson 
line so that it is minimized instead by configurations in 
which tr [Q n ] = 0, which in turn implies that the eigen- 
values of £1 are uniformly distributed around the unit 
circle. Thus, the set of eigenvalues is invariant under the 
Ztv transformations, which multiply each eigenvalue by 
e 2mk/N Rotate the unit circle by k/N). The center sym- 
metry is then unbroken by construction. The coefficients, 
{a n }, can be suitably chosen such that the deformation 
potential, P[fi], forces unbroken symmetry at any com- 
pact ificat ion scale [17], but for our purposes we are only 
interested in small compactifications (L <C A -1 where L 
is the length of the compactified dimension and A is the 
QCD scale). At small compactification, the gauge cou- 
pling at the compactification scale is small so that the 
semiclassical computations are under complete theoreti- 
cal control [17]. 
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B. Infrared description 

As discussed in [17], the proper infrared description (at 
distances larger than the compactification scale) of the 
theory is a dilute gas of N types of monopoles, charac- 
terized by their magnetic charges, which are proportional 
to the simple roots and affine root of the Lie algebra for 
the gauge group U(1) N . The extended root system is 
given by the simple roots, 



(5) 



ot-i = (1,-1,0, ...,0) = ei -e 2 , 

a 2 = (0, 1, -1, ... ,0) = e 2 - e 3 , 

&N-1 = (0, . . . ,0, 1, -1) = e N -i-e N , 
and the affine root, 

a N = (-1,0, . . . ,0, 1) = ejsr-ei. 



We denote this root system by A a ^ and note that the 
roots obey the inner product relation 



a a • a b = 2S a ,b - £ a ,6+i - S a . 



6-1- 



(6) 



For a fundamental monopole with magnetic charge 
a a G A a ff, the topological charge is given by 



Q 



-J 



1 



16tt 2 



tr 



= ± 



TV' 



(7) 



and the Yang-Mills action is given by 



Sym 



•'xS 1 



1 

2? 



tr [F%] 



(8) 



/ 

Jr 



d 4 x 



R3XS 1 



2g 2 



rtr 



8tt 2 



IQI- 



The second equivalence hold because the classical 
monopole solutions are self dual, F^ v = ±F^ U . 

The ^-parameter in the Yang-Mills action can be in- 
cluded in conventional way, 



Sym —> Sym + i0 



with = e^ pa F ni 



1 



167T 2 



tr 



(9) 



The system of interacting monopoles, including 6 pa- 
rameter, can be represented in the dual sine-Gordon form 
as follows [17, 22], 



$dual 



C 



(10) 



where ( is magnetic monopole fugacity which can be ex- 
plicitly computed in this model using the conventional 
semiclassical approximation. The parameter enters 
the effective Lagrangian (10) as 0/N which is the di- 
rect consequence of the fractional topological charges of 



the monopoles (7). Nevertheless, the theory is still 2ir 
periodic, but not because of an explicit 2tt periodicity 
of Lagrangian (10). Rather, it is restored as a result of 
summation over all branches of the theory when the levels 
cross at 6 = ir (mod 2tt) and one branch replaces another 
and becomes the lowest energy state as discussed in [22] . 

In the following sections we shall need an explicit ex- 
pression for the topological density and magnetic field in 
terms of scalar a field, 



S(x) 



16tt 2 



tr 



F^F^ 



-1 



N 



Ajk4 ' 



N 



4tt 2 



,(a) 
M 



V-B^(x) 



(a) „(a) 
jk r iA 



(ii) 



where the U(1) N magnetic field, B l = e tjM Fj k /2g is ex- 
pressed in terms of scalar magnetic potential as follows 



2ttL 



e ijk d k *W, 



2ttL 



Va (a) . (12) 



The expression for the magnetic field in terms of scalar 
magnetic potential should not be surprising as our system 
is in fact magnetostatic and a description in terms of cr^ 
is quite appropriate to study the relevant dynamics. 

The explicit form for the creation operator for a 
monopole of type a at x is given by [22] 



r(x) 



(13) 



A^a(x) = J 

and for an antimonopole by 

M a (x) = e - iaa,<T(x) . (14) 

The expectation values of these operators (A4 a (x)) in 
fact determine the ground state of the theory. Formula 
(13) shows again that <x(x) can be interpreted as a mag- 
netic scalar potential. 

Finally, the dimensional parameter which governs the 
dynamics of the problem is defined as 



(15) 



This parameter can be interpreted as Debye correlation 
length of the monopole's gas. The average number of 
monopoles in a "Debye volume" is given by 



AT--" 3 



m„ 



<-(£)' 



1 



»i, 



(16) 



The above inequality holds since the monopole fugacity 
is exponentially suppressed, ( ~ e -1 ^ , and we should 
view (16) as a constraint on the validity of the approxi- 
mation where semiclassical approximation is justified. 
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III. DOMAIN WALLS IN DEFORMED QCD 

There is a discrete set of degenerate vacuum states as 
a result of the 2tt periodicity of the effective Lagrangian 
(10) for the a field, and thus there exist domain wall 
configurations interpolating between these states. The 
corresponding configurations are not however conven- 
tional domain walls similar to the well known ferromag- 
netic domain walls in condensed matter physics which 
interpolate between physically distinct vacuum states. 
Here, instead, the corresponding configuration interpo- 
lates between topologically different but physically equiv- 
alent winding states |n), which are connected to each 
other by large gauge transformation operator. There- 
fore, the corresponding domain wall configurations in Eu- 
clidean space are interpreted as configurations describing 
tunnelling processes in Minkowski space, similar to Eu- 
clidean monopoles which also interpolate between topo- 
logically different, but physically identical states. 

In fact, a similar domain wall which has an analogous 
interpretation is known to exist in QCD at large temper- 
ature in weak coupling regime where it can be described 
in terms of classical equation of motion. These are so- 
called Zn domain walls which separate domains charac- 
terized by a different value for the Polyakov loop at high 
temperature. As is known, see the review papers [27, 28] 
and references therein, these Zn domain walls interpolate 
between topologically different but physically identical 
states connected by large gauge transformations similar 
to our case. These objects can be described in terms of 
classical equation of motion and have finite tension ~ T 3 
such that their contribution to path integral is strongly 
suppressed. While the corresponding topological sectors 
are still present in the system at low temperature (though 
they are realized in a different way) it is not known how 
to describe the fate of Z^r walls within QFT in the strong 
coupling regime wherein the semiclassical approximation 
breaks down. 

The domain walls to be discussed below in deformed 
QCD are very much the same as Zn domain walls at 
high temperature and their contribution to path integral 
is also strongly suppressed as their tension is finite in 
weak coupling regime. Nevertheless, one can study the 
structure of these domain walls, as well as their inter- 
action with dynamical magnetic monopoles. We con- 
jecture that the domain walls we describe below in the 
weak coupling regime in deformed QCD slowly become 
the objects (with effectively vanishing tension) which are 
observed in lattice simulations [1-8] in the strong cou- 
pling regime, as we adiabatically increase the coupling 
constant without hitting the phase transition as argued 
in [17]. This portion of the theory can not be tested in 
our deformed QCD model in the semiclassical approx- 
imation, but hopefully this portion of strongly coupled 
dynamics can be understood in the future using different 
techniques, such as the dual holographic description as 
advocated in the present context in [13, 14]. 



A. Domain wall solution 

There are many different types of domain walls sup- 
ported by the system (10) which have very different phys- 
ical meaning. In this paper we focus on the discrete sym- 
metry of the effective Lagrangian (10) given by the 2tt 
shift, a a — > a a + 27r, where any component of a field 
can be shifted by 2ir independently. There are N dif- 
ferent domain wall types similar to monopole case since 
classification of our system is based on cti G A a ff . We em- 
phasize that there are only (N — 1) physical propagating 
photons in the system as one scalar singlet field, though 
it remains massless, completely decouples from the sys- 
tem, and does not interact with other components at all 
[17]. As a result of this structure, a configuration with 
N different types of magnetic monopoles will carry zero 
magnetic charge and one unit of the topological charge 
Q = 1 as each monopole carries Q = 1/N topological 
charge. The corresponding configuration can be identi- 
fied with a conventional instanton with Q = 1 which is 
made of N constituents. A similar comment also applies 
to DW structure: a configuration with N different types 
of DWs on top of each other will produce a trivial vac- 
uum configuration as the (N — 1) abelian components of 
the magnetic field will cancel each other, similar to mag- 
netic monopole construction. Thus, although there are 
N different types of the DWs in our construction, only 
(TV — 1) of them are independent. 

In what follows, without loosing any generalities, we 
consider N = 2 case. In this case there is only one phys- 
ical field x which corresponds to a single diagonal com- 
ponent from the original SU (2) gauge group. The action 
(10) becomes, 




In terms of x field, the classical equation of motion which 
follows from (17) and which determines the profile of the 
domain wall has the form, 

V 2 x-m^sinx = 0, (18) 

where we take = to simplify things, and the mass of x 
field m x = 2m a is related to the Debye correlation length 
(15). The solution of this sine-Gordon equation which 
interpolates between x( z = —oo) = and x( z = +oo) = 
27r, and which is centered at zo = being independent of 
x, y coordinates is well known 

x(z) = 4arctan [exp(m x z)] . (19) 

As we mentioned before, the domain wall (19) does not 
describe a physical domain wall (DW) which interpolates 
between physically distinct vacuum states, but rather in- 
terpolates between topologically different but physically 
identical states. We remark that a similar construction 
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Topological Charge Density 



FIG. 1. This picture explains the transition between paths 
corresponding to the decay of some domain wall state to a 
domain wall free ground state. The path wrapping the peg 
represents a state with some domain walls, while the path 
that does not denotes a state with no domain walls We can 
deform the DW path by lifting it over the obstacle so that 
we can unwind it and deform it into the DW-free path. If 
the path describes domain walls with some weight, then it 
would require some energy to lift over the obstacle. If this en- 
ergy is not available, then classically, the configurations that 
wind around the peg are stable. Quantum mechanically, how- 
ever, the domain wall could still tunnel through the peg, and 
so the configurations are unstable quantum mechanically, see 
estimate for this probability in Appendix A. Picture adapted 
from [31]. 



has been considered previously in relation with the so- 
called N = 1 axion model [29, 30], more recently in the 
QCD context in [31], and in high density QCD in [32]. 
In the previously considered cases [29-32] as well as in 
present case (19) there is a single physical unique vacuum 
state, and interpolation (19) corresponds to the transi- 
tion from one to the same physical state. Therefore, such 
domain walls are not stable objects, but will decay quan- 
tum mechanically, see Appendix A for corresponding es- 
timates. Nevertheless, if life time of configuration (19) 
is sufficiently large, it can be treated as stable classical 
background, and it can be used to study the interaction 
of domain walls with monopoles, which is one of the main 
objectives of present work, see Figure 1 with more expla- 
nations. 

The width of the domain wall is determined by m^ 1 ' 
while the domain wall tension a for profile (19) can be 
computed and it is given by 



+oo 



L 2 



(0 



L 3 ' 



(20) 



With explicit solution at hand (19), the magnetic field 
(12) distribution inside the domain wall is given by 



B-, 



= (JL-\ 



4m y 



V47rL/ (e r - 



(21) 




FIG. 2. This picture shows the two layer structure of the 
topological charge density. 



The topological charge density distribution can then be 
computed using formula (11) with the following result 



q(z) - -sm X - 77 ( e m x * + e -m x *)2 



(22) 



From equation (22), we see that the net topological 
charge Q ~ dzq(z) on the domain wall vanishes. 
However, the charge density has an interesting distribu- 
tion; it is organized in a double layer structure, which 
is precisely what apparently has been measured in the 
lattice simulations [1-4]. For a graphical depiction see 
Figure 2. The same double layer structure can be seen 
by computing the magnetic charge density pM which is 
defined as 



(a) 

Pm 



v - B(a) «M^) 



dz 2 



<) 2 



(23) 



Thus, the relation between the topological charge den- 
sity (22) and magnetic charge density (23) holds for the 
domain wall 



\LN J 



(24) 



in full agreement with the general expression (11). 

From eqs. (22), (23), we see that an average density of 
magnetic monopoles filling the interior of domain wall is 
expressed in terms of the same parameter £ which char- 
acterizes the average monopole's density in the system 
(10). One can interpret this relation as a hint that the 
topological charge sources have a tendency to reside in 
vicinity of the domain walls rather than being uniformly 
distributed. We further elaborate on this matter in sec- 
tion IV. 

It is interesting to note that the domain walls in de- 
formed QCD model are very similar (algebraically) to the 
well known domain wall studied previously in some SUSY 
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models, see e.g. review [33]. Of course, there are funda- 
mental differences between the two: in SUSY models the 
domain walls interpolate between physically distinct vac- 
uum states, in huge contrast with our domain walls which 
correspond to interpolation between topologically differ- 
ent but physically identical states. Therefore, the inter- 
pretation in these two cases is fundamentally different: 
in SUSY models the domain walls are real physical ob- 
jects, while in our deformed QCD model they should be 
interpreted similar to instantons, objects which describe 
the tunnelling processes, see [13] with more comments 
on this interpretation. Furthermore, the classification of 
the domain walls in SUSY models is based on flavour 
group symmetry breaking SU(Nf) — >> [/(l)^ -1 , in con- 
trast with colour symmetry breaking in deformed QCD. 
However, the formal classification of the domain walls in 
SUSY models based on simple roots from flavour group is 
very much the same as classification in our case based on 
SU(N) [/(l)^- 1 breaking pattern, see (5), (6). These 
similarities include, in particular, highly nontrivial prop- 
erties such as ordering of the domain walls or their pass- 
ing through each other. However, these questions will 
not be elaborated on in the present work. 

The most important lesson from this analysis is that 
the double layer structure naturally emerges in the con- 
struction of the domain walls in the weak coupling regime 
in deformed QCD. As claimed in [17] the transition from 
high temperature weak coupling regime to low temper- 
ature strong coupling regime should be smooth without 
any phase transitions on the way. Therefore, it would be 
tempting to identify the double layer structure found in 
this work (22) with the double layer structure from lat- 
tice measurements [1-4] when one slowly moves along a 
smooth path from the weak coupling to the strong cou- 
pling regime. 



IV. DW-MONOPOLES 

We now consider the domain wall configurations dis- 
cussed in the previous section interacting with monopole 
configurations. In these computations the domain walls 
are treated as classical background fields, and as such we 
do not consider fundamentally quantum questions, such 
as the density of domain walls. Instead, we consider some 
questions which can be answered in the semiclassical con- 
text. We focus on the interaction between a monopole 
and domain wall, each acting as magnetic sources, and 
compute the energy of the configuration as a function 
of separation distance between the two. Therefore, the 
question we are addressing is where would a point charge 
prefer to sit in the presence of our domain wall? We 
should our domain walls are not empty objects as they 
are already filled by magnetic monopoles with density 
determined by (22). 

Again, we consider the simplified scenario of S77(2), 
which corresponds to considering the interaction between 
a single type of domain wall, a, and a monopole of the 



same type, or to be more precise an antimonopole so that 
the magnetic charge is — a a . The domain wall is defined 
as previously, (19), but centered at a distance zq from 
the origin, so that the magnetic scalar potential is given 
by (letting m = m x ) 



Xz ( x ) = 4 arctan 



o m(z-z ) 



(25) 



The monopole is defined such that it is a point source 
solution to the Klien-Gordon equation, 



vV 



mtp = 0, 



(26) 



centered at the origin, and is thus an approximate solu- 
tion to the sine-Gordon, (18), away from the origin. The 
magnetic potential of the monopole is then given by the 
well known Yukawa potential, 



<p(x) = 



-m|x| 



rax 



(27) 



We then consider the configuration of monopole and 
domain wall separated by a distance zo and would like to 
compute the magnetostatic energy (Euclidean action) as 
a function of z$. The energy associated with just a do- 
main wall alone is proportional to the area of the domain 
wall, which is infinite in this case, so we compute instead 
the difference between the energy of the two together and 
the energy of the two independently, 

AE (z ) = S [ Xzo +<p]-S [ Xzo ] -S[<p], (28) 

where S is given by (axes have been rescaled relative to 
(17)) 



six] 



/ d 3 x 

JR3 



\ (Vx) 2 - rn 2 cosx 



(29) 



The quantity AE is finite but we cannot compute it an- 
alytically, and so we compute above integrals numeri- 
cally instead, for zo varying near the domain wall. Some 
technical details of the computation are as follows. We 
work in a cylindrical volume around the origin with ra- 
dius 10/ra and length 30/ra. It is valid to do so since 
the monopole potential is exponentially suppressed with 
length constant ra. We were forced to remove a small 
volume around the origin when computing the potential 
energy term because the structure is that of the sine of 
a divergent quantity, which is highly oscillatory. These 
two approximations make up the bulk of the uncertainty, 
which is ~ 5%, for most values of Zq. When comput- 
ing the values of zq with the monopole center near a 
peak in the charge density of the domain wall, the inte- 
gral converges much more slowly, and so we turned down 
the working precision, which in turn increased the uncer- 
tainty to ~ 20%. This was however only required for five 
points near the peak and trough of the charge density of 
the domain wall, that is zq « 1/ra (see Figure 2). 

Performing the numerical integration results in the plot 
given in Figure 3. There is an attractive potential be- 
tween the monopole and domain wall with the monopole 
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Energy as a Fuction of Seperation Distance 



m 




Displacement 



FIG. 3. This is a plot of the numerical result of change 
in energy computed at various separation distances between 
domain wall and monopole. Notice that for zo < 0, the 
monopole to the right of the domain wall, there is an "at- 
tractive" potential with a minimum at ~ l/3m. 



on one side (zq < 0), and a slightly repulsive one for 
the other side (zo > 0). Also, there is a minimum at 
zq ~ l/3m , so that the monopole would prefer to sit in- 
side the domain wall, between the center and the peak of 
the sheet with the same charge density. It is both inter- 
esting and unexpected that the monopole (with charge 
—a) is attracted to the domain wall sheet with the same 
charge (—a) rather than the sheet of opposing charge 

(«)■ 

Figure 3 is not the complete story since we have not 
considered possible changes in the magnetic flux distribu- 
tion coming from the monopole. Basically, the monopole 
shape could deform in response to the interaction with 
the domain wall, so as to become less spherically symmet- 
ric. In order to properly treat this problem, we should 
allow the spherical distribution of the monopole to vary 
to some superposition of solutions to the Klien-Gordon 
equation (26). This described further calculation is be- 
yond the scope of this work, but we conjecture that the 
magnetic field will prefer to orient itself along the do- 
main wall, so that the magnetic flux will be pushed out 
to the edge of the domain wall at the boundary of space, 
similar to arguments presented in refs. [19-21] in cosmo- 
logical context. In this way, we have a picture in which 
any point-like magnetic monopoles become bound to ex- 
tended domain walls with any magnetic flux being pushed 
along the domain walls to infinity. Apparently, this is 
precisely the picture discovered in lattice simulations [1- 
4] when very few localized 4d objects are observed in the 
system, see item 6) in the Introduction. 



V. CONCLUSION, SPECULATIONS, AND 
FUTURE DIRECTIONS 

There are two important results of this work. Firstly, 
a double layer structure similar to that which is observed 
in lattice simulations [1-4] naturally emerges in the con- 



struction of the domain walls in weak coupling regime in 
deformed QCD. Secondly, monopole configurations char- 
acterized by well localized topological (and magnetic) 
charge interact with domain walls in such a way that 
there is an attraction between the two, and the monopole 
favors a position inside the domain wall. We introduced 
these domain walls as external background fields, while 
they are expected to be dynamical configurations with 
effectively vanishing tension in strong coupling regime as 
holographic picture suggests. We further believe that the 
magnetic field due a monopole in the presence of a do- 
main will tend to align with the domain wall, such that 
the flux is pushed to the boundary of the domain wall. 
If this effect persists in strongly coupled regime, it could 
be an explanation for the observation in lattice simula- 
tions [1-4] that there are no well localized objects with 
finite size which would carry the topological charge. As 
to weather the described realignment of magnetic field in 
the presence of a domain wall is energetically favoured is 
an interesting future problem. 

In weak coupling the domain wall solution is a nicely 
behaved smooth function, but what happens when we 
transition slowly to the strong coupling regime? The 
holographic picture suggests that the effective domain 
wall tension vanishes and so they can be formed easily 
in vacuum. It is possible that the domain walls become 
"clumpy" with a large number of folders. Such fluctu- 
ations would then increase the entropy of the domain 
wall, which eventually could overcome the intrinsic ten- 
sion. If this happens, the domain walls would look like 
very crumpled and wrinkled objects with large number of 
foldings, and as such, the domain walls may loose their 
natural dimensionality, and become characterized by a 
Hausdorff dimension as recent lattice simulations suggest 
[9]. Nevertheless, the topological charge distribution on 
larger scales after averaging over a large number of these 
foldings should be sufficiently smooth so that the double 
layer structure would not disappear because the transi- 
tion from weak to strong coupling regime should be suf- 
ficiently smooth as argued in [17]. Therefore, we identify 
the double layer structure found in this work (22) with 
the double layer structure from the lattice measurements 
[1-4]. These particularities of the transition from weak 
to strong coupling are also interesting future questions, 
which will likely require an analysis beyond the semi- 
classical level. 

Finally, it would be interesting to see how the effects 
which we studied in weakly coupled deformed QCD man- 
ifest themselves in the dual holographic description [14]. 
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Appendix A: Decay of domain walls 

The decay mechanism is due to a tunnelling process 
which creates a hole in the domain wall which connects 
the x = domain on one side of the wall to the x — 2tt 
domain on the other, see (19). Because the ground 
state on the two sides is physically identical, it is pos- 
sible for the fields to remain in the ground state as they 
pass through the hole. That is, there is no interpolation 
(winding) as they pass through the hole. This lowers the 
energy of the configuration over that where the hole was 
filled by the domain wall transition by an amount propor- 
tional to R 2 where R is the radius of the hole. The hole, 
however, must be surrounded by a string- like field con- 
figuration which interpolates between an unwound con- 
figuration and a wound one. This string represents an 
excitation in the heavy degrees of freedom and thus costs 
energy, however, this energy scales linearly as R. Thus, 
if a large enough hole can form, it will be stable and 
the hole will expand and consume the wall. This process 
is commonly called quantum nucleation and is similar to 
the decay of a metastable wall bounded by strings; there- 
fore, we use a similar technique to estimate the tunnelling 
probability. The idea of the calculation was suggested in 
[29] to estimate the decay rate in the so-called N = 1 
axion model. In QCD context similar estimations have 
been discussed for the rf domain wall in large TV QCD in 
[31] and for the rf domain wall in high density QCD in 
[32]. 

If the radius of the nucleating hole is much greater 
than the wall thickness, we can use the thin-string and 
thin- wall approximation. This approximation justified as 
we shall see when we calculate the critical radius R c . In 
this case, the action for the string and for the wall are 
proportional to the corresponding worldsheet areas 



S (R 3 x 8 1 ) = 2tt RLa - irR 2 Lcr. 



(Al) 



The first term is the energy cost of forming a string, 
where a is the string tension and 2ttRL is its worldsheet 
area. The second term is energy gain by the hole over 
the domain wall, in which a is the wall tension and ttR 2 L 
is its worldsheet volume. We should note that formula 
(Al) replaces following, more familiar expression for the 
classical action which was used in many previous similar 
computations, see [31, 32] 



S (R 4 ) 



4ttR 2 c 



4-7T 

~3~ 



(A2) 



Minimizing (Al) with respect to R we find the critical 
radius R c and the action So 



R c = -, S K 3 x S 1 = , 

a v ' a 



(A3) 

which replace more familiar expressions for the critical 
radius R c = ^ and classical action Sq(M 4 ) = 1Q %JZ from 
[31,32]. 

Therefore, the semiclassical probability of this process 
is proportional to 



T ~ exp 



ira 2 L 



(A4) 



where a is the DW tension determined by (20) , while a is 
the tension of the vortex line in the limit when the inter- 
action term ~ ( due to the monopole's interaction in low 
energy description (10) is neglected and U(l) symmetry 
is restored. In this case the vortex line is a global string 
with logarithmically divergent tension 



a ~ 2ir 



— (-Y 

AL 2 \2ttJ 



In 



R 



R c 



(A5) 



where R ~ m" 1 is a long-distance cutoff which is deter- 
mined by the width of the domain wall, while R cor e ~ L 
when low energy description breaks down. The vortex 
tension is dominated by the region outside the core, so 
our estimates for computing a to the logarithmic accu- 
racy are justified. Furthermore, the critical radius can 
be estimated as 



Rc — 



2m, 



ln( 



1 



(A6) 



which shows that the nucleating hole ~ R c is marginally 
greater than the wall thickness ~ m" 1 as logarithmic fac- 
tor ln(^-j-) ~ In A/" ^> 1 where Af ^> 1 is large parame- 
ter of the model, see (16). Therefore, our thin-string and 
thin- wall approximation is marginally justified. 

As a result of our estimates (A4), (20), (A5) the fi- 
nal expression for the decay rate of the domain wall is 
proportional to 



T ~ exp 



7ra 2 L 



exp — 7T 



3 ( JM y ™x LJ 



exp (-7 -Af In 2 Af) <C 1, 



(A7) 



with 7 being some numerical coefficient. The estimate 
(A7) supports our claim that in deformed QCD model 
when weak coupling regime is enforced and Af ^> 1 the 
domain walls are stable objects and our treatment of the 
domain walls as stable objects is justified. 
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